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ABSTRACT 

Observations with the Infrared Space Observatory (ISO) SWS spectrometer are used to con- 
strain the velocity law and wind clumping of the well-studied Wolf-Rayet (WR) star WR 136 
(HD 192163) (WN6). Because the free-free continuum opacity in WR winds increases steadily 
with wavelength in the IR, each point in the continuous spectrum may be regarded as forming in 
a pseudo-photosphere of larger radius for longer wavelength. Using this idea in combination with 
an analysis of the Doppler-broadened widths of several He II recombination lines, we can derive 
information about the velocity law and dumpiness of the stellar wind of WR 136. The observed 
line emission emerges from the region exterior to the radius of optical depth unity in the free-free 
opacity, corresponding to v > O.Svoo for our shortest wavelength line. The ISO observations 
provide the continuum shape, flux level, and seven fairly strong He II emission profiles. Adopting 
a /3-law distribution for the outflow velocity law, we compute the continuous energy distribution 
and line profiles. We find that there is a broad range of /3-values consistent with the continuum 
data if we also allow the wind temperature to be a free parameter. Interestingly, the continuum 
data are found to constrain the wind to have fairly low clumping values for the IR-forming re- 
gion. Using the continuum results in conjunction with line profile modeling, the observational 
constraints are best satisfied with a clumping factor of Dc —< > / < p >^ oi 1-3 and a /3- value 
of 2-3, although higher /3 values are not strongly ruled out for a modest wind temperature. The 
wavelength range of our ISO data allows us to probe only the outer wind acceleration zone, but in 
combination with radio observations, our finding that the wind clumping is fairly small suggests 
that the clumping in the wind of WR 136 decreases with increasing radius. 
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1. INTRODUCTION 

The Wolf-Rayet (WR) class of early-type stars 
has long been known for having extreme stellar 
wind properties. In the visible band, broad emis- 
sion lines that form in the wind dominate the stel- 
lar spectrum (Abbott & Conti 1987). The expla- 
nation of these extreme outfiows has long posed se- 
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vere problems for stellar wind theories. In partic- 
ular these stars are noted for their "wind momen- 
tum problem" (Barlow et al. 1981; Abbott et al. 
1986; Willis 1991), in which the momentum flux 
in the wind Mvoo (where M is the mass loss rate 
and Voc is the terminal speed of the stellar wind) 
can exceed by more than an order of magnitude 
the momentum flux L^/c (where is the stellar 
Bolometric luminosity) of the radiation field pur- 
portedly responsible for the outflows. This situa- 
tion is normally parametrized by the performance 
number 77 — Mv^oc/ L^, which is of order 10 for 
WR stars (although there can be a wide range - 
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sec Hamann & Koesterke 1998). In recent years 
the observational estimates of WR wind momenta 
have been reduced, partly from upward revisions 
of stellar luminosities (by line blanketing, e.g., 
Hillier & Miller 1999) and partly by accounting for 
dumpiness in the winds (Hillier 1991; Hamann & 
Koesterke 1998), yet values of Mvoo remain large 
(« 10, Nugis & Lamers 2000). 

Alternative wind-driving mechanisms have 
been proposed to explain the large wind mo- 
menta, but several studies have indicated that the 
WR winds can be driven by radiative forces alone, 
without invoking other physics (such as magnetic 
fields), through multi-line scattering of the stellar 
radiation field (Lucy & Abbott 1993; Springmann 
1994; Gayley, Owocki, & Cranmer 1995; Onifer 
& Gayley 2003). Multi-line scattering is possible 
because of ionization gradients in the wind flow. 
As photons escape from the inner wind with an 
opacity characterized by a particular distribution 
of lines, the photons encounter different sets of 
line opacities at larger radii because of changes in 
the dominant ions. Indeed a wide range in ioniza- 
tion states as a function of radius appears to exist 
in the WR winds (Kuhi 1973; Schultc-Ladbeck, 
Eenens, & Davis 1995; Herald et al. 2000). Thus 
there can be a greater deposition of momentum by 
the stellar radiation field to drive the wind, and 
an enhancement of the wind momentum flux. 

Although multi-line scattering can influence the 
mass-loss rates of some O stars, it tends to be 
much less important than for the WR stars (Vink, 
de Koter, & Lamers 2000). As we argue below, 
one of the consequences of multi-line scattering in 
the wind is that the rate at which the wind ve- 
locity increases with radius should be shallower 
in WR winds than in typical O stars. Observa- 
tional constraints on the velocity laws of WR stars 
should thus provide a test of the multi-line scat- 
tering wind acceleration process. Obtaining such 
information has been the motivation for this ob- 
servational study. 

For early type stars, it is traditional to describe 
the wind velocity distribution as a function of ra- 
dius using the "/3 velocity law", which typically 
has the form, 

vir)=v^(l-^Y . (1) 



This formula for the velocity distribution was de- 
rived with /3 = 0.5 in one of the earliest of 
radiation-driven outflow models by Milne (1926). 
The same value was predicted in the line-driven 
wind theory of Castor, Abbott, & Klein (1975; 
CAK). In detailed modeling of Of and OB super- 
giants by Grocnewcgen & Lamers (1989), values 
of P =0.7-0.9 were found. These higher values in- 
dicate a shallower rise in the velocity distribution 
versus radius and are also in better agreement with 
the predictions of the modified CAK (mCAK) the- 
ory of Friend & Abbott (1986) and Pauldrach, 
Puis, & Kudritzki (1986). These mCAK models 
accounted for the finite size of the star produc- 
ing the wind. They also led to predictions for M 
and Vaa that are in good agreement with obser- 
vations of O stars. The difference in the /? val- 
ues in going from the CAK to the mCAK model 
arises from the fact that the angular distribution 
of the driving radiation field is different. For the 
very forward-peaked radiation field of CAK, the 
acceleration is rapid, but for the somewhat more 
dispersed radiation field associated with the fi- 
nite disks of mCAK, the velocity law is shallower. 
The multi-line scattering loads to an even less 
strongly forward-peaked radiation field and to yet 
another step away from a point-like radiation field. 
Thus it is to be expected that the velocity law for 
WR stars should have even higher values of (3. 

With the view that multi-line scattering can 
lead to a significant enhancement of the winds, 
momentum transfer is no longer the major is- 
sue in understanding the large mass-loss rates of 
WR stars. Instead there is now a need for the 
opacity to adequately "trap" or prevent the leak- 
age of the photons so that multi-line scattering 
can indeed occur. As the photons transmit some 
of their outward momentum to the wind, they are 
red shifted. If there are gaps in the opacity versus 
wavelength distribution, the photons can escape 
and thus contribute no further to the wind ac- 
celeration. Since the winds of WR stars require 
that photons transmit about 77 x their momentum 
hvo /c to the wind, a very large effective number of 
scatterings per photon is required. For example, 
Gayley et al. (1995) show in their Figure 3, that a 
photon will interact with approximately t/^ lines, 
and experience a net redshift of AA » Aq {r]Voo/c). 
A photon that initially interacts with a EUV line 
near 600 A will tend to drift about 66OA redward. 
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Consequently, there can be few gaps in the hne 
opacity distribution, lest photons flood through 
opacity minima to escape the wind before deposit- 
ing sufficient momentum to the flow. Hence the 
wind momentum problem has now come to be in- 
terpreted as an "opacity problem" (Gayley et al. 
1995). 

Rather little is known observationally about 
the radial distribution of velocity in WR winds. 
However, a good estimate for the wind terminal 
speed Woo can be made from forbidden line pro- 
files (e.g.. Barlow, Roche, & Aitken 1988; Willis 
et al. 1997; Morris et al. 2000; Ignace et al. 2001 
[hereafter Paper I]). This is because such lines 
are typically formed at large distances of hun- 
dreds of stellar radii or more, where the wind has 
achieved its asymptotic expansion speed. Because 
the forbidden lines are optically thin and form in 
the constant velocity zone, they tend to c;xliibit 
flat-topped profiles, the half- widths of which ac- 
curately determine the wind terminal speed. 

Much less is known about the radial distribu- 
tion of the velocity v{r) deeper in the wind. There 
have been several discussions about the values of 
(3 or the consequences of shallow velocity laws 
in WR stellar winds. From theoretical modeling, 
Schmutz (1997) arrived at a, (3 ranging from 3 to 
8. Monte Carlo radiation transport simulations 
by Springmann (1994) and Gayley et al. (1995) 
also reveal extended wind acceleration zones (not 
characterizable by a single /? value). Support for 
these theoretical results comes from the study of 
binary systems by (Auer & Koenigsberger 1994). 
Also, Lepine & Moffat (1999) have derived large f3 
values from observations of the regular motion of 
small emission peaks or ''migrating features" that 
appear in WR He II emission lines. An interest- 
ing structural effect of large /? values was disciissed 
by Ignace, Cassinelli, & Bjorkman (1996). They 
found that velocity laws with modestly large val- 
ues of /3 (such as /3 = 3) are more likely to form 
wind-compressed zones as a result of stellar rota- 
tion. This effect was used to explain why some 
WR stars show intrinsic polarization (e.g., Har- 
ries, Hillier, & Howarth 1998). 

To study the velocity law in WR stars further, 
we obtained Infrared Space Observatory (ISO) 
satellite observations, and in this paper we de- 
scribe results from the best data set, which cor- 
responds to WR 136. We use observations of both 



the IR free-free continuum and emission profiles of 
He II recombination lines. In §2, we discuss in gen- 
eral the continuum and line diagnostics of the den- 
sity and velocity distributions in WR stars. In §3, 
the description of the data obtained with the ISO 
satellite and its analysis to infer /3 is presented. 
The results of the study are discussed in §4. 

2. INFRARED DIAGNOSTICS FOR DENSE 
STELLAR WINDS 

Relative to the flows from other early-type 
stars, the WR winds are denser and they produce 
spectra that are dominated by emission lines. The 
understanding of these spectra involves three key 
elements: (a) In the optical, a major continuum 
opacity is electron scattering, which is gray but 
involves complicated thermalization effects (Miha- 
las 1978) (b) In the UV and EUV spectral ranges, 
strong line blanketing both blocks the flow of ra- 
diation and drives the wind, (c) In the IR through 
to the radio spectral region, it is the increasingly 
large free-free opacity that accounts for the IR and 
radio continuum excesses. The advantage of using 
this long wavelength spectral region to study the 
wind structure is that the continuum opacity in- 
creases monotonically with wavelength. For wind 
temperatures of a few times 10^ K that are appli- 
cable to radii where the free- free continuum forms, 
the opacity varies as in the near IR, gradu- 
ally shifting to a dependence at longer wave- 
lengths. In the latter case, the free-free opacity 
is so great that the emergent radiation originates 
at large radii where the flow has reached terminal 
constant expansion-velocity. The free-free opac- 
ity is lower in the IR, and this permits radiation 
to emerge from the parts of the wind where the 
flow is still accelerating. In particular, ISO ob- 
servations provides us with the spectral data that 
allows us to diagnose an interesting segment of 
the wind velocity law. There are two complemen- 
tary approaches that we employ: (a) a study based 
on the continuous energy distribution, and (b) a 
study of emission line profiles formed by recombi- 
nation. The latter method is, in fact, related to 
the first in that the continuos opacity of the wind 
affects the range of depths over which we receive 
line emission, and thus can have a strong influence 
on profile strengths and shapes. Consequently, we 
consider the continuum first, and then proceed to 
the emission line profiles. Many of the issues rel- 
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evant to the interpretation of IR observations of 
WR stars have been discussed in investigations 
similar to ours, such as Hilher, Jones, & Hyland 
(1983) and HiUicr (1987ab) in detailed studies of 
WR 6 (= HD 50896 = EZ CMa). 

2.1. IR Continuum Slope Methods 

The classic approach to studying WR density 
and velocity distributions is to consider the con- 
tinuum slope. Wright & Barlow (1975) and Pana- 
gia & Fclli (1975) wore the first to demonstrate 
that the spectral slope should vary as oc 
(or fx oc A~^'^) for a spherical wind in constant 
expansion-velocity and thick in the free-free con- 
tinuum. There are several informative scalings 
that came from these early studies and subsequent 
considerations. For example, Wright & Barlow 
showed that in the constant expansion-velocity 
case, the continuum emission is equivalent to a 
volume integral over the cmissivity oc with a 
lower bound Rv,yo\ determined by the radius at 
which the line-of-sight free-free optical depth is 
Tjy.ff = 0.244. This optical depth is somewhat 
smaller than the 2/3 value associated with the 
Eddington-Barbier result for plane-parallel stel- 
lar atmospheres. Adopting an "equivalent photo- 
sphere" approach, Hillier et al. (1983) showed that 
for the observed monochromatic luminosity with 
= 47ri?^g^j.f7ri?,y(ro), where is the Planck 
function at the wind temperature, the monochro- 
matic radius Rv,smf corresponds to a radial depth 
of T^.ff = 0.05. This small optical depth means 
that to understand the monochromatic luminosity, 
it is essential to account for the spatially extended 
nature of the wind continuum formation region. 

As in Wright & Barlow (1975) and Hillier et al. 
(1983), let us compare the line-of-sight optical 
depths for different wind velocity laws. For the 
geometry depicted in Figure 1, the free- free opti- 
cal depth along a ray that passes through the en- 
velope at impact parameter p (normalized to i?*) 
is given by 



ni{x)ne{x) dz' , (2) 



for X = ^ z'^ + = r/Rt: the normalized stellar 
radius, z = Z/R^ the line-of-sight coordinate to- 
ward the observer, m the ion number density, 



the electron number density, and (T, v) a fac- 
tor that isolates the temperature and frequency 
dependence of the free-free opacity. This opacity 
factor is given in Wright & Barlow as 



Jfff (T, v) = 3.7 X 10^ g^{T) T'^^ y-^ 

X (l - e-^"'^'^) cm^ Hz(^) 

where Z{ is the rms ion charge, gu{T) is the Gaunt 
factor, T is in K. and v is in Hz. We assume that T 
is constant throughout the continuum-formation 
portion of the wind, and its value is treated as a 
free parameter, for which we know some limiting 
values. 

We assume a spherical radial wind and use the 
dimensionless radius x = r/R and velocity w = 
v/voo- From the continuity equation, the radial 
distribution of density is given by 



M 



Po 



47r r^u(r) x^w{x)' 
where the scale factor is 



M 



Po 



(4) 



(5) 



and the normalized wind velocity distribution is 



w{x) = 11 — 



/3 



for 



b = 



1 - 



(6) 



(7) 



Equation (6) differs from equation (1) in that we 
allow for the parameter 6 < 1 and provides for 
a finite initial value wq = vq/voc at the base of 
the wind where x = 1. To evaluate the free- free 
optical depth in equation(2), the electron and ion 
densities are related to the mass density as rie = 
p/ ideUiH and Ui = pj iJimH, for /Xe and /xi the mean 
molecular weight per free electron and free ion, 
respectively. 

For illustrative purposes, let us compare two 
simple cases: constant expansion- velocity (i.e., 
/? = 0) and a /? = 1 velocity law involving typical 
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wind acceleration versus radius. The comparison 
is motivated by the fact that we seek to use the 
observed spectral energy distribution to infer in- 
formation about the wind velocity law. With these 
two values of /3 we gain insight into how the veloc- 
ity law influences the observables, and conversely 
what our ISO observables can tell us about the 
velocity distribution. 

For the constant expansion-velocity case, the 
line-of-sight optical depth Tf[{x) = Tf[{z,p = 0) 
becomes 



from an assumed wind velocity law toward a con- 
tinuum slope, the relation between the wind ex- 
pansion and observed continuum emission can be 
addressed with an "inverse approach" . Cassinelli 
& Hartmann (1977) considered the relation be- 
tween the observed spectral slope and the density 
and temperature variations assuming power law 
forms, with p oc and T oc r"™. Using the 
effective radius concept, they derived a relation 
between the power-law slope vfi, = Xfx oc A* and 
the q and m indices as given by 



(8) 



so that a specified optical depth is obtained at 
a radius x-r = (tc/3t)^/^, where the constant 
collects all of the various wind, atomic, and fre- 
quency dependencies in a single parameter. The 
total radial optical depth down to the wind base 
at .T = 1, is Tc/3. The emergent intensity from any 
point along the line-of-sight is exp(— Tj^^ff ). 

For the (3 = 1 velocity law, the integral for the 
optical depth is again analytic, giving 



T-c(i^) 
63 



1 



1 - bu 



- {1 - bu) + 2ln{l - bu) 
(9) 

where u = 1/x. The total line-of-sight optical 
depth to the wind base is Tc b~^[b{2 — b)/{l — b) + 
21n(l -6)]. 

In Figure 2, both the optical depth and the 
emergent intensity Ii, are plotted against radius 
for these two cases. In this example both Kg and 
po have been held constant. There are four im- 
portant points to note, (a) The total line-of-sight 
optical depth is greater for larger /3, because at 
every radius, the flow will be slower and hence 
denser, (b) A fixed optical depth occurs at larger 
radius for larger (3. (c) Thus the observed free-free 
emission emerges from larger radius for larger /?, 
leading to a greater monochromatic free-free lumi- 
nosity L^. And (d) for a constant source function 
Bu{T), intensities emerge from contributions over 
a broad range of radii. When evaluating the to- 
tal flux from the unresolved source, the net effect 
of larger values in (3 is to steepen the wind den- 
sity distribution, which in turn leads to a steeper 
continuum spectrum. 

Instead of the "forward problem" of working 



s = 



6q 



7 



-2q 



1 



(10) 



Taking the wind to be isothermal with m = 0, 
equation (10) reduces to s = {6q — 7)/(l — 2q). 
For a spherical wind, inferring q for the density 
law is equivalent to determining the shape of the 
wind velocity distribution, although not the val- 
ues of the velocity. Furthermore, given the range 
in wavelengths of the IR observations, the shape 
of the derived velocity distribution is determined 
over the corresponding range in continuum forma- 
tion radii. More generally, the spectral slope will 
depend not only on the density distribution, but 
also on the temperature gradient, electron scatter- 
ing, bound-free opacity, and wind ionization vari- 
ations. The main point of this discussion is that 
the wind velocity law can control the continuum 
spectral slope and may thus be inferable from ob- 
servations. 

2.2. The IR Emission Line Profile Method 

Our data set contains a number of IR recombi- 
nation lines of He II. The strong lines observed in 
our spectra are for transitions involving levels in 
the range 6 to 16. Assuming pure recombination 
lines and accounting for free- free opacity, we model 
the line profiles using standard Sobolev theory in 
spherical symmetry (e.g., Mihalas 1978). The to- 
tal emergent intensity arriving at the Earth along 
a ray at p (see Fig. 1) is 



(11) 



where the first term on the RHS is the emergent in- 
tensity of the underlying continuum emission and 

the second term is the emergent intensity of line 
emission. Both the line and continuum emission 
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arc computed at each frequency in the Une profile. 
In Sobolev theory, the fine emission at a given fre- 
quency arises from an "isovelocity zone" , which is 
a region of constant Doppler shift in the flow as 
perceived by the observer. Figure 1 shows an ex- 
ample of an isovelocity surface. The continuum 
emission, however, is determined by a volume in- 
tegral over the whole envelope. The continuous 
opacity can attenuate the line emission; moreover, 
the line opacity can diminish the net continuum 
emission. 

It is convenient to define four optical depths. 
The first three arc for the continuos opacity: Tforc, 
Taft, and Ttot- Of these the first is the optical 
depth that attenuates the line emission (i.e., by 
matter that lies on the near side of the isoveloc- 
ity zone with respect to the observer). The sec- 
ond does not (i.e., matter on the rear side of the 
isovelocity zone). The third is simply the sum, 
Ttot = Tfore + ^aft- The fourth optical depth is the 
Sobolev optical depth rg associated with the line 
opacity. This is given by 



TS = 



\dvz/dz\ 



4 2/ \ 
X W [X) — — 

dz 



(12) 



where the line-of-sight velocity gradient is 



2 dw 

dx 



= M- — + (13) 



dWr, 

dz ax ' X 

where fj, = cos 9, and the scale constant ti is 



Ki po X i?* 



(14) 



Using these optical depth definitions, the emer- 
gent continuum and line intensities are given re- 
spectively by 

/■not /"Tforc 

/.,c= / B^T) e-^ e-^^ dT+ B^{T)e^^dT 

(15) 

and 



(16) 



If we assume that the wind is isothermal, or at 
least approximately constant over the radii giving 
rise to the bulk of the line and continuum emission, 
the intensity for the continuum reduces to 



(17) 

for Tw the wind temperature. If the further as- 
sumption is made that the departure coefficients 
for the upper level populations are nearly unity, 
one has that S^.i = Bu{T^), and the total emer- 
gent intensity for line and continuum emission, in- 
cluding the relevant absorptions, along a ray at p 
compactly reduces to 



h{p) = B,{T^ 



1-e 



-(not+Ts) 



(18) 



Clearly, outside the line frequencies, the calcula- 
tion for the continuum emission is given by the 
same expression only with = 0. We ignore any 
bound- free contribution to the opacity (sec the dis- 
cussion of §4) , but note that the bound- free emis- 
sivity at IR wavelengths also scales as so that 
the form of equation (18) would continue to hold. 
Electron scattering is also ignored, but it becomes 
more important at shorter wavelengths where the 
free-free opacity is smaller. Equation (18) implies 
that in the line, the intensity is always less than 
that in the neighboring continuum. In such a case 
the only way that a line will appear in emission 
is if the integration of the intensity over the ap- 
parent disk of the envelope is larger at line wave- 
lengths. Thus an emission line profile provides us 
with information about the velocity distribution 
somewhat beyond the radii at which the neigh- 
boring continuum is formed. 

We are interested in using the line data to learn 
about the wind velocity distribution. Figure 3 
shows the emergent flux of line emission, given by 
ji, exp(— T^^fj) dV, for the line emissivity and 
dV a differential volume element, along the line- 
of-sight toward the center of the star. The upper 
panel is a plot against continuum optical depth, 
and the lower is against the normalized wind ve- 
locity. The two solid curves at top are for explic- 
itly optically thin lines, contrasting the constant 
expansion-velocity case with /? = against the 
(3 = 1 case. Also plotted is a dotted line for an 
optically thick line with /3 = 1. From Figure 3, we 
may conclude that the bulk of the line emission 
emerges exterior to the radius where ri,,ff = 1. In 
the bottom panel, only the thin (solid) and thick 
(dotted) line cases for /? = 1 are shown, since there 
is no variation of v with r for /? = 0. We see that 
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the line emission originates from a fairly broad 
range of physical wind velocities. At longer wave- 
lengths, where the free- free opacity is larger, the 
peak of the flux contribution will move to higher 
wind velocities resulting in emission lines of in- 
creasing breadth. 

To obtain the total line flux at a given line fre- 
quency (or velocity shift v = v^) requires an inte- 
gration over impact parameter for the correspond- 
ing isovelocity zone, as given by 

/. = ^ / ^-^P^P'^P' (19) 

where d is the distance to the star. For thin 

lines the integrand of equation (19) scales with 
/O^r^ exp(— Ty^ff) oc x^^w"^ exp(— Ty^ff). The width 
of the line profile depends on the wind expansion 
in the region where the line forms. Since there is 
a finite range of radii that primarily contribute to 
the line emission, the overall width of the line is set 
by the velocity at the radius where Tff « 1. Line 
emission from smaller radii is strongly absorbed, 
and that from larger radii is small owing to the 
diminishing emission measure. 

In this discussion of the line radiative trans- 
fer, the forward approach has been emphasized. 
As was the case for the continuum, it is possi- 
ble to use an inverse approach to analyze emission 
lines, if certain simplifying assumptions hold. The 
inverse approach can allow one to directly deter- 
mine the wind velocity law (Cannon 1974; Brown 
et al. 1997; Ignace et al. 1998ab). The primary as- 
sumptions arc that the lines be optically thin, the 
wind be spherically symmetric, and that there be 
a relatively distinct photospheric radius. Unfortu- 
nately, consideration of the S/N, difficulties with 
line blends, non-thermal broadening (e.g., from 
wind instabilities), and the radial extension of the 
free-free continuum formation zone, have led us 
to conclude that the application of these inverse 
techniques to the observed He II profile shapes in 
WR 136 will not be useful. The method of Ignace 
et al. (1998b) appears to show some promise. It 
uses the distribution of total line intensities for 
several thin lines, and the total emission for each 
line depends on the emission measure. Thus some 
of the effects that severely bias the inversion of line 
profile shapes are averaged out. However, there 
are not enough suitable lines in our dataset to at- 



tempt this inversion. 

Returning to the forward approach, a key point 
in modeling the IR spectrum of dense winds is that 
both the line and continuum optical depths scale 
with the square of the density. The assumption 
that the source function for the line and contin- 
uum are the same ensures that there is no net 
absorption feature in the profile. So the situa- 
tion results in rather symmetric emission profiles, 
which initially is somewhat counterintuitive, since 
the absorption of line emission by the continues 
opacity varies strongly from the near side of the 
star (corresponding to blueshiftcd line frequencies) 
to the far side (for redshifted line frequencies). 

A counter example is provided by wind X-ray 
emission lines that are expected to be asymmet- 
ric in the standard paradigm of distributed wind 
shocks (Ignace 2001; Owocki & Cohen 2001; Ig- 
nace & Gayley 2002). For X-ray lines, the emis- 
sivity scales with p^, and the line emission suffers 
photo-absorption by the dominant "cool" wind 
component. Asymmetry in the profile shape re- 
sults because the continuous opacity makes no 
contribution to the emission, and such asymmetric 
morphologies are observed in the Chandra spectra 
of C Pup (Cassinelli et al. 2001). Indeed in the case 
at hand, we have examined the line emission con- 
tribution alone (i.e., with attenuation by free-free 
absorption) and find it to be strongly asymmetric 
with blueshiftcd peak emission, just like the theo- 
retical X-ray lines. However, the absorption of the 
continuum is also asymmetric, with greatest ab- 
sorption occurring at blueshifted frequencies. This 
is why the net profile recovers a fairly symmetric 
appearance. 

To show profile effects explicitly, an analytic ex- 
pression for the profile shape can be derived in 
the case of a constant expansion-velocity flow, if 
the free-free optical depth is quite large. In fact, 
the profile shape in this limit has previously been 
derived by Hillier et al. (1983), but the result ap- 
pears there as an unnumbered equation (between 
their eqs. [13] and [14]), and examples of the profile 
shapes are not shown. Those authors were more 
concerned with line equivalent widths. (In their 
appendix a derivation for the equivalent width is 
presented for a wind in homologous expansion, 
with V (X r.) The solution for this special case 
is particularly apt to the problem at hand, so we 
provide details for the line profile derivation in the 
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appendix. Results are shown in Figure 4, where 
a sequence of hne profiles are plotted as the ratio 
of line to continuum optical depth scales, ti/tc, 
is varied. The profile is rather "bubble-shaped" 
and is identically symmetric. Of course, relax- 
ing any of our assumptions will, in general, allow 
the profiles to develop asymmetries; however, the 
bubble-shape profiles illustrated by the constant 
expansion- velocity case are similiar to those actu- 
ally observed in WR 136, as we describe next. 

3. ISO OBSERVATIONS OF WR 136 

We have obtained ISO observations of the pro- 
totype WN6 Wolf-Rayet star WR 136 (HD 192163). 
As a member of the Cyg OBI association, the dis- 
tance to this star is well-determined [d— 1, 820 pc; 
Lundstrom & Stenholm 1984). WR 136 is unusual 
among WN6 stars in that it has a substantial 
abundance of hydrogen, ^ 12% by mass (Crowther 
& Smith 1996; Hamann et al. 1994). Line and 
continuum spectra of WR 136 were obtained with 
the SWS spectrometers aboard the ISO satellite 
as described in Table 1. Automated preliminary 
processing of the data was done using version 7 
of the standard pipeline package (OLP-7), and 
subsequcmt processing was done as described in 
detail in Paper I. In this section we present the 
continuum slopes and line profile widths that can 
be used to derive density and velocity information 
for the wind of the star. 

3.1. The Continuum Flux Distribution 

Continuum measurements were made over the 
range of 2.8 to 10 microns, beyond which the S/N 
ratio was insufficient to reliably determine a con- 
tinuum. The spectrum obtained with the SWSOl 
spectrometer is displayed in Figure 5, and this has 
been dereddened using the algorithm of Whittet 
(1992), with values for Ay and Eb-v given by 
van der Hucht (2001). A power-law fit was made 
to the continuous flux distribution, as shown in 
Figure 6, that yields a power-law index of 



If 



d log fx 
dlogX 



-2.87. 



(20) 



This value is to be compared with a slope of —2.6 
expected from a constant expansion wind at long 
wavelengths (Wright & Barlow 1975). The steeper 



observed value suggests that cither the wind is ac- 
celerating, the temperature is varying, or both. 
However, we find that the power law is somewhat 
weakly constrained, especially by the low quality 
data at longer wavelengths. The la-confidence 
interval for acceptable power-law indices is quite 
asymmetric, with a range of —2.75 to —2.9 that 
applies to the better quality short wavelength 
spectrum (the —2.87 is the minimum value). 
This range is consistent with spectral indices re- 
ported by Nugis, Crowther, & Willis (1998), who 
quote a value of — 2.75±0.02 in the range of 12 /^m 
to 6 cm, and a value of —3.15 ± 0.35 in the range 
of 12- 25 ^m. 

Using the range of power laws from our ISO 
data, we can estimate the density variation over 
the range of radii in which the continuum forms 
for these wavelengths. In the power-law approach 
of Cassinelli & Hartmann (1977) described in 
§2.1, the s in equation (10) is related to 7ir 
via 8 — 7f + 1, which in our case ranges from 
— 1.75 to —1.9. Assuming an isothermal wind with 
m = 0, the density power law is found to be 
q={s^ 7) /{2s -h 6), which ranges from 2.1 to 2.3. 
But of course, if p oc r""?, and p oc r~^w~^, then 

V oc r*~^, which in this case results in u w r"'^. 
This is a fairly shallow slope, and consequently 
the continuum data, based on this crude analysis, 
is likely arising from rather far out in the wind 
acceleration zone, where the wind is starting to 
approach constant expansion- velocity. 

In terms of a standard wind [3 law like equa- 
tion (6), one can make a further identification. 
Assuming large radius with r ^ R^,, the velocity 
gradient is d\nv/d\nr « (ibR^Jr. With v oc r°'^, 
we have that (3 « 0.2r/6i?*, which implies that 

V « Woo(l — /36i?*/r) w 0.8woo regardless of the 
/3-cxponent. This exercise suggests that our data 
probe mainly the outer wind acceleration. 

3.2. The He II Line Profiles 

The high resolution SWS06 observing mode was 
used to obtain recombination line data with re- 
solving powers ranging from A/AA » 1300 to 2500 
over wavelengths of 2.8 to 5.0 microns, respec- 
tively. We use the seven best He II recombination 
lines in our study. Figure 7 shows the six lines used 
for determining the line half-widths (the seventh, 
14-11 is too noisy to obtain a reliable width, but 
is adequate for determining an equivalent width). 
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Line identifications with comments arc given in 
Table 2. Atomic data and measurements for these 
lines are provided in Table 3. The HWHMs for 
the hnes are around 1000 kms"-'^, and the instru- 
mental broadening is about 100 kms^^. Since the 
emission profiles are somewhat Gaussian, estimat- 
ing the effect of instrumental "smearing" is rel- 
atively straightforward. The convolution of two 
Gaussians yields another Gaussian. The width of 
the resultant Gaussian is given by the rms of the 
widths for the two input Gaussians, in this case 
the emission line and the instrumental smearing. 
Thus, instrumental broadening affects our HWHM 
determinations at about the 1% level, and so we 
have ignored it. 

Using the line width data, we can make an- 
other estimate of the wind velocity law in WR 136. 
HWHM values were measured for the red wing of 
the six He II recombination lines and are plotted 
against wavelength in Figure 8. (HWHMs mea- 
sured from the blue wing are systematically less 
by about 100 kms~^, which wc attribute to the 
influence of line absorption.) A power-law fit was 
made to this data, yielding a power-law index of 



7v = 



d l og HWHM 
dlogX 



= 0.22 ±0.07, (21) 



To infer information about v{r), a power-law in- 
dex slope in radius is needed. For this, one could 
assume that the line emission originates from the 
vicinity of the pseudo-photosphere of the contin- 
uum (i.e., around Ti,fi = 1), since the line emission 
interior to this region will be strongly absorbed. 
But without knowing the wind velocity law a pri- 
ori, it is impossible to compute the radius of opti- 
cal depth unity in the free-free opacity. We can at 
least estimate the wavelength dependence of the 
radius as follows. 

From the observed continuum emission, we can 
define the surface equivalent radius that would 
provide the observed flux, with 



Rx 



surf 



= d 



fx 



(22) 



Already from §2.1, we have argued that this ra- 
dius i?A,surf will be larger than the radius where 
Tiy^ff = 1. Even so, the overall wavelength de- 
pendence of -RA.surf should roughly mimic that of 
the radius of optical depth unity, namely that at 



longer wavelengths, the two radii will increase in 
size in a manner that is roughly proportional. 

With fx and Bx taken as known, the run 
of i?A,surf with wavelength can be determined, 
and we find that = dlogRx/dlogX = 2 + 
0.5 J F or about 0.6. So we may infer that 
dlogw/(ilogi?A,surf = Iv/jR ~ 0.4, which is con- 
sistent with the estimate (of 0.2) derived from 
§3.1 based purely on the continuum slope data. 
Of course, these two approaches (continuum slope 
and line width) are not entirely independent, since 
7i{ is related to jf- Nevertheless, the line- width 
data is independent, and so the consistent results 
support the argument that we are observing the 
outer wind acceleration at these wavelengths. 

3.3. Combined Analysis of Continuum and 
Lines 

We have written a code to compute both free- 
free emission and recombination line profiles using 
the formalism of §2.2. We use this code to repro- 
duce the continuum fluxes and line profile shapes 
for our best He II lines. The goal is to simultane- 
ously constrain the velocity law /?, the characteris- 
tic wind temperature T^, and the wind clumping. 

3.3.1. Clumping Factor Effects 

Clumpiness, now a well-recognized attribute of 
many WR stars (e.g., Moffat et al. 1988; Hillier 
1991; Nugis, Crowther, & WiUis 1998; Hamann & 
Koesterke 1998), is included in our modeling. We 
define the clumping factor parameter as 



(23) 



The wind clumping is assumed to occur on small 
scales, so as to produce the same emission mea- 
sure as a smooth wind but with a smaller M. 
The clumping may occur on a range of scales, but 
for our theoretical approach, these scales must be 
small enough that the Sobolev approximation is 
not invalidated. In essence, this means that the ra- 
diative transfer may still be described by integral 
expressions, in contrast to discrete sums (Hamann 
& Koesterke 1998b). 

3.3.2. Assumed Wind Parameters 

We consider the terminal speed to be known 
from Paper I, at v^o = 1490 kms~^. For the mass- 
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loss rate, Tabic 4 lists several values for WR 136 
as taken from the literature; the method of evalu- 
ation; and whether the estimates assumed smooth 
or clumpy winds. The first four values assumed 
an underlying smooth wind. For the last value, 
Nugis et al. (1998) used radio observations to de- 
rive M values using a clumped wind model, and 
found that at the large radii where the radio emis- 
sion is produced, the wind of WR 136 is consistent 
with no clumping (i.e., Dc ~ 1). On the other 
hand, Crowther & Smith (1996) conclude in fa- 
vor of clumping in the wind of WR 136 because 
otherwise their models overestimate the electron 
scattering wings of strong optical/NIR lines. That 
these two different methods yield different clump- 
ing predictions must reflect on the fact that the 
methods sample different radial locations in the 
wind. It may not be surprising that clumping can 
be a fimction of radius (e.g., as in the theoretical 
models of Runacres & Owocki 2002), and at least 
for WR 136, the degree of clumping seems to di- 
minish in going from the inner wind to the outer 
wind. 

It is not our goal in this study to determine the 
mass-loss rate; instead, the unclumped mass-loss 
rate from Nugis et al. is adopted. To study the 
influence of clumping for the IR continuum and 
line formation, we assume that the clumping factor 
is constant over the radii in which these emission 
form. Our model fits to the observations can thus 
place constraints on Dc- 

Also needed for our models is a value of the 
"stellar radius", i?*. For other classes of stars, 
this would typically be the radius of the wind base. 
For example in the case of Of and OB stars, it is 
possible to observe a photospheric radiation field 
from which one can derive basic quantities associ- 
ated with a "hydrostatic atmosphere" , such as the 
stellar rotation speed, surface gravity, and effec- 
tive temperature. Consequently, there is a well- 
defined radius associated with the wavelength at 
which a hot star exhibits peak emission (i.e., in 
the UV/EUV). The only free parameters in the 
canonical /3 velocity law equation (1) are the value 
of P itself and the terminal velocity Woo • The lat- 
ter is generally determined from the strong and 
saturated P-Cygni lines in the UV (e.g., Prinja, 
Barlow, & Howarth 1990) or the widths of forbid- 
den lines. 

However, the situation is different for WR stars, 



fundamentally because the base of the wind is 
not observed at any wavelength (Abbott & Conti 
1987). The mass-loss rates are so high for the 
WR stars that the winds arc optically thick at 
all wavelengths. This gives rise to the well-known 
problem that for such stars, the meaning of stel- 
lar radius and effective temperature arc ambigu- 
ous (e.g.. Castor 1970; Cassinelh 1971). There 
have been two approaches to finding values to use 
for i?* and Tofr in WR stars. One is that model 
atmospheres are computed that have thick winds 
(e.g., Hamann & Koesterke 1998ab). The stel- 
lar radius is chosen to be at the depth at which 
the Rosseland mean optical depth reaches a value 
of 20. Using this radius and the stellar lumi- 
nosity, a stellar temperature is defined from 
the Stefan-Boltzmann relation. The values for 
WR 136 derived in this way are = 70800 K 
and i?* = 7.5i?0 (Hamann & Koesterke 1998a). 

Nugis & Lamers (2000) use another approach 
based on the idea that there must be a definable 
hydrostatic star at the base of the wind. These 
authors start with an estimate of the radius iJevoi 
of WR stars provided by stellar evolution mod- 
eling. For example, using formulae derived by 
Schaerer & Maeder (1992) for H-deficient stars, 
Nugis & Lamers find for WR 136, i?evoi = 1-2Rq 
and Tevoi = 140600 K. Their radius is 1/6*'* as 
large and their temperature 2 times larger than 
values derived with the wind modeling approach. 
Since WR 136 is not completely H-deficient, the 
radius derived from the evolution models is prob- 
ably underestimated; however, it should still be 
smaller than that from the wind-modeling method 
above. 

While the i?cvoi approach may provide a good 
estimate of the underlying size of the star, there 
is a problem in that the optical depth through the 
wind to the central star is extremely large, orders 
of magnitude greater than unity. Given that the 
radiation field is not strongly peaked in the out- 
ward direction deep in the wind, we expect that 
the rate of increase of the outflow speed is likely to 
be shallow just above i?evoi- Hence, we do not con- 
sider it likely that this evolutionary radius is the 
proper one to use in a /3-law distribution that de- 
scribes the conditions in the IR line and continuum 
formation regions of the wind. In other words a ve- 
locity law prescription governed by a single value 
for f} is inadequate for describing both the slowly 
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accelerating inner wind and that in the outer wind. 
Nevertheless, it is convenient to use the /3-law for- 
mula to discuss the outer part of the wind that 
wc can probe with our data. Thus the b R prod- 
uct that appears in the velocity law, would not 
in this case correspond to the evolutionary radius 
of the star, nor would the velocity, vq necessarily 
be the sound speed at that radius, but are simply 
parameters that complete the description of the 
outer velocity law. 

To bring closure to our models, we choose to 
adopt the value i?* = 5Rq from the wind mod- 
eling of WR 136 by Crowther & Smith (1996). 
Their model is also useful for limiting the pos- 
sible range in wind temperature. For WR 136, 
Crowther & Smith (1996) find an effective tem- 
perature of 27700 K at a Rosseland optical depth 
of 2/3, but a "stellar" temperature of 57700 K at 
an optical depth of 20 corresponding to a stellar 
radius of i?* = 5Rq. In modeling the forbidden 
lines at very large radius (lOO's of i?*), an electron 
temperature of 13000 K is used (e.g., in Paper I). 
Recent models of WR winds show variations in 
with radius, with being large at the in- 
ner wind, then decreasing with radius to achieve 
an asymptotic value (e.g., Dessart et al. 2000 and 
Herald et al. 2001). In this paper we are attempt- 
ing to extract information from the ISO observa- 
tions about the velocity distribution without doing 
full radiatively-driven hydrodynamic simulations 
that might ultimately be required. It is likely that 
the insight developed from this empirical analy- 
sis will provide guidance to more detailed studies. 
So as in Hillier (1987ab), we take the wind to be 
isothermal in the IR formation region of interest 
to us. We take the temperatures 27700 K and 
57700 K as reasonable limiting values in our mod- 
els. The basic stellar and wind parameters used 
for WR 136 are summarized in Table 5. 

3.3.3. The Fit to the Continuum 

For the continuum emission in the vicinity of 
our line measurements, we have computed a grid 
of solutions which are shown in Figure 9. There 
are two sets of plots: one for = 27700 K and 
one for 57700 K. The three different curves shown 
as solid in the figure are for /3 = 1,3, and 10 (from 
left to right). The different points are for different 
clumping factors, with Dc = 1,3,10,20, and 50, 
such that Dc = 1 is the lowest point on each curve. 



The horizontal and vertical dotted lines demarcate 
"allowed zones" from the observations, represent- 
ing the IfT confidence range of models permitted 
by the errors in the measurements. For the power- 
law fits to the continuum slope, we use the form 
fx = /o {^/^oY, with Ao = 1 micron. As quoted 
previously, values of p range from —2.75 to about 
—2.9, with —2.87 giving the minimum fit- The 
absolute flux error for ISO is about 15%, and prop- 
agating errors for the variation in p results in the 
allowed range for the flux scale. Favored models 
are those that match both the observed continuum 
slope (plotted as the abscissa) and the overall flux 
scale (plotted as the ordinate), and fall within the 
indicated error box. Note that we do not actually 
observe the spectrum at 1 /im; rather the 1 /im 
flux level /o is a back-extrapolated fiducial value 
based on the power-law fit to the spectrum. 

For comparison we also show models in which 
the He ionization changes from He III to He II at 
a radius of r = 30i?* (the dashed line curves). We 
chose this radius for the change in ionization on 
the basis of model calculations by Hillier (1987a) 
for the well studied WN 4 star WR 6. There 
are probably differences that will arise by the fact 
that WR 136 is a WN 6 star instead of WN 4, 
and because of the fact that WR 136 is not hy- 
drogen free. However, we want to investigate the 
possible effects of the ionization change and find 
that Hillier's study is the best currently available 
guide. The location of the ionization transition re- 
quires the detailed models to provide, for example, 
the ionizing EUV flux that penetrates through the 
wind. We suspect that the radius of the ionization 
transition in WR 136 must be relatively large, be- 
cause we observe such strong He II recombination 
lines. 

The effect of an ionization change is to shift 
the curves of Figure 9 primarily in a lateral fash- 
ion, and more so for the lower temperature models 
than for the higher temperature ones. These new 
curves would alter conclusions as to the velocity 
law. However we note that rather low clumping 
factors are still required in the fit shown in the 
figure. Clearly, low clumping factors are found 
both with and without the ionization transition. 
To investigate the other effects that infiuence the 
velocity law diagnostics, we shall assume for the 
remainder of the paper that the wind is fully dou- 
bly ionized He. 
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Focusing on the solid lines of Figure 9, the 
observed continuum flux data favor a /3-value of 
about 5 for the cooler wind, and 3 for the hotter 
wind. These values are larger than the value /? = 1 
commonly assumed in WR wind models. However, 
values as low as /? = 1 are not strongly ruled out 
in Figure 9. The influence of the temperature can 
be understood as follows. A lower temperature re- 
duces the free-free source function, but increases 
the free-free opacity. For a wind that is isother- 
mal and in constant expansion, as applies to the 
radio emission, the temperature dependence ex- 
actly cancels out (CassincUi & MacGregor 1986). 
This cancellation is no longer the case in the ac- 
celerating portion of the wind. Also, the effect 
of the wind acceleration is to modify the density 
distribution. In the net a larger (3 both elevates 
and steepens the density distribution. Increasing 
(3 thus leads to a brighter continuum level with a 
steeper power-law slope. 

Consider the effect of changing the clumping 
factor. An inc;rease in the clumping factor raises 
the flux level. However for a fixed (3, the enhanced 
clumping elevates the density such that the ob- 
served continuum emission forms at larger radius 
and higher speed flow. Consequently, the contin- 
uum slope is actually reduced, since the continuum 
forms in a region closer to constant expansion. 
Thus to match both the fiducial flux level and the 
continuum slope, larger /3-values are needed for 
winds of lower temperature. 

Note that in the lower temperature case, the 
power-law slope can actually fall below the canon- 
ical 0.6 value of Wright & Barlow (1975) for a 
constant expansion- velocity wind. This occurs for 
two reasons; (a) both and A are sufficiently 
low that the Planck function is not entirely in the 
Rayleigh-Jeans limit, which influences the slope, 
and (b) the wavelength dependence of the Gaunt 
factor in the IR band is slightly different from that 
in the radio band. For the IR free-free Gaunt fac- 
tor, we use data from Waters & Lamers (1984). 
These authors computed Gaunt factors down to 
10/j.m only; however, we adopt their expressions 
for the shorter wavelengths of interest with He^"*" 
as the dominant contribution (with He mass frac- 
tion Y = 8/9). We also include free-bound opac- 
ity from H+ and He^^, again using the Waters 
& Lamers expressions. The free-bound opacity 
from the IR and longward has the same scaling 



as the free-free, and so the effective Gaimt factor 
becomes g^, = gs + bib, where 6fb is generally much 
smaller than gg- 

3.3.4- Fits to the Line Profile Data 

To further constrain wind parameters, we have 
also produced model line profiles to be compared 
with the observations. To model the lines with the 
influence of wind acceleration, we again assume 
that He++ is the dominant ionization state and 
calculate the populations of the upper He+ levels 
by assuming LTE with the continuum. The NLTE 
departure coefficients 6k for the upper levels of in- 
terest (fc > 7) are expected to be near unity (e.g.. 
Storey & Hummer 1995). A grid of line profile cal- 
culations are shown in Figure 10, with variations 
in /3, line optical depth t;, and the wind tempera- 
ture Tvy. Although these parameters are not truly 
independent - for example, r; will depend on both 
/? and T„, they are treated as independent for the 
purpose of indicating the range of profile shapes 
that can result. The "canonical" model to be used 
as a reference in this figure (shown as dashed) is for 
/3 = 1, Dc = 1, Tw = 30000 K, and n = 1. Varia- 
tion of Dc is not explicitly shown, but its influence 
is similar to temperature, in the sense that increas- 
ing the clumping or decreasing the wind temper- 
ature both result in overall weaker emission lines 
relative to the continuum level. Note that in our 
line modeling, we take the inner wind density to be 
fixed, with /9(i2*) = M/AttE^vq a constant. This 
means that as (3 is varied, h also changes to keep 
vq fixed, with vq = 0.0075t;oo- 

Before applying the line fits to the data, we 
first demonstrate that the observed lines are sam- 
pling part of the wind acceleration by overplot- 
ting the asymptotic expression for emission line 
profiles with the observations in Figure 11 (see 
Eq. [A15]). For our seven best lines, the emis- 
sion profile shapes from a wind expanding at the 
terminal speed are consistently too broad, and so 
the line formation must occur in the lower velocity, 
accelerating portion of the flow. 

In fitting the observed profiles, both (3 and 
are allowed to vary; however, we must also al- 
low for additional broadening to approximate the 
influence of "turbulence" in the wind, and so a 
"broadening parameter" a is introduced. This 
additonal broadening is motivated by the com- 
monly held fact that the winds of hot stars are 
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structured. The velocity distributions of hot star 
winds are non-monotonic, with variations associ- 
ated with wind shocks that reach outflow speeds 
that can be significantly greater than the asymp- 
totic value obtained at large radius. Line-driven 
wind instability effects can produce sharp drops 
and spikes in the wind velocity and density dis- 
tribution in O star winds (Lucy & Solomon 1970; 
Lucy 1982; Owocki, Castor, & Rybicki 1988; Feld- 
mcicr et al. 1997). Similar effects are theoreti- 
cally expected for WR winds (Gayley & Owocki 
1995), and structure is inferred from UV and opti- 
cal emission profile variations (e.g., St-Louis et al. 
1995; Lepine & Moffat 1999). Of specific rele- 
vance to this paper, an analysis of the C IV dou- 
blets near 1550 A reveal blueshifted black troughs 
and somewhat larger absorption edge velocities for 
many classes of hot star winds, including WR stars 
(Prinja et al. 1990). The difference between the 
trough and edge velocities, amounting to as much 
as 25% of Voo for WR winds, is commonly inter- 
preted as arising from the wind being structured 
in velocity. For WR 136, Prinja et al. report an 
edge velocity of 2700 kms~^ but a terminal speed 
of 1600 kms~^, amounting to a difference of over 
1000 kms~^. Based on our previous comment 
about the convolution of Gaussian profile shapes, 
broadening of this order would increase the line 
widths of our model profiles by about 40%. 

To constrain values of /? and for the IR 
lines, we have computed a grid of line pro- 
files to fit our best line. He II 7-6, which is 
strong and has excellent S/N. Using the termi- 
nal speed of Woo = 1490 kms~^ from Paper I, we 
compute the required additional broadening from 
a = sjHWHMl^^ - HWHM'^^^, with results 
shown in Figure 12, to match the line HW (defined 
as the line HW in the red wing to be compared 
with the observed values from Tab. 3). The solid 
lines arc for the hotter wind, and the dashed ones 
for the cooler wind temperature. In each case, 
the upper curve is for = I and the lower for 
Dc = 3. Larger values of Dc at a fixed (3 will 
produce broader line widths resulting in lower cr- 
values, but the analysis of the continuum emission 
(based on Fig. 9) rules out Dc greater than about 
3. We did consider the influence of a transition in 
the He-ionization at r = 30i?*. With this transi- 
tion the continuum- normalized line emission drops 
by about 10% or less, and the shape is not dra- 



matically affected. The effect is not major because 
both the line and free-free opacities are affected 
so as to be partially compensating when the line 
is continuum normalized, and so we continue to 
assume that the wind is pure He III. 

In Figure 12, the lower dotted line is the in- 
strumental broadening, representing a minimum 
threshold for a. A strong upper limit would be 
1200 kms~^ based on the UV C IV line absorp- 
tion (not shown in the Figure). Here we have an 
even better constraint from the fact that the wind 
cannot have clumping factors below unity. At this 
point, the lines would suggest that /3- values below 
1 are probably not consistent with the observed 
7-6 fine. 

However, we can do better than this. Wolf- 

Rayet stars are sources of X-ray emission (e.g., 
Pollock 1987; Wessolowski 1996; Ignace, Oskinova, 
& FouUon 2000). Recent pointed observations of 
3 putatively single WR stars by XMM-Newton 
indicate hot gas of at least a few million Kelvin 
in the winds (Skinner 2002ab; Ignace, Oskinova, 
& Brown 2003). If arising from wind shocks, as 
would be the standard explanation, shock velocity 
jumps of at least a few hundred kms~^ would be 
needed to generate hot gas at these temperatures. 
A range of 300-500 kms~^ would seem to be rea- 
sonable, and such variations relative to the mean 
flow would lead to "broadening" of the IR lines. 
The observed X-rays do not emerge from the same 
radii as the IR continuum and He II lines; the X- 
rays emerge mainly from farther out in the wind, 
but only because of strong photo-absorption by 
the dense WR wind. 

Consequently, we indicate in Figure 12 with 
a circled asterisk a "best representation" for the 
wind of WR 136 that takes into account (a) low 
clumping based on the continuum analysis, (b) the 
rather broad but finite range of /3 values from the 
continuum constraints, and (c) a modest level of 
a necessary to account for the fact that WN stars 
are X-ray sources. The result is = 2, with 
/? = 3 for the cool wind model and /? = 2 for the 
hot wind model. Naturally, based on the discus- 
sion of the He ionization, and ambiguity regard- 
ing the broadening processes intrinsic to the winds 
of WR stars, (3 is somewhat loosely constrained, 
although a clumping factor that is not too large 
seems required. 

Using these preferred values for /? and Fig- 
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ures 13 and 14 show the model profiles (dotted) as 
overplotted on the observed emission lines (solid), 
for the cool and hot wind cases respectively. Note 
that the lines in these figures arc continuum nor- 
malized. A key for the line transitions is pro- 
vided at bottom center. In addition to broadening 
processes, line blends may influence some of the 
observed profiles widths, particularly in the wide 
He II 16-10 line and the wings of He II 9-7, and 
the somewhat blueward broadened He II 7 6 line. 

Inspection of our model results for the lines also 
confirms the earlier estimates that the ISO data 
are not sampling the innermost wind acceleration. 
Optical depth unity in the free-free opacity occurs 
at a radius of > O.Swoo for our shortest wave- 
length line. Still, our analysis places interesting 
constraints on the value of /3 and for the ex- 
tended wind acceleration of WR 136. 

4. DISCUSSION 

We have argued that WR stars should have 
a significantly different velocity distribution than 
OB stars, for which the radiation field is more 
forward-peaked near the base of the wind. The 
many thick lines and the small mean free paths 
deep in the star should cause the initial rise in 
velocity to be gradual. One would expect that a 
flow which is driven by multi-line scattering would 
have a large value of /3, a conclusion that we have 
sought to test using IR observations of WR 136 
obtained with the ISO satellite. Line and contin- 
uum data have been analyzed and are found to be 
sensitive primarily to the outer part of the wind 
at the observed wavelengths. Although our ISO 
data do not probe the deep layers where multi- 
line scattering processes are most important, the 
velocity range that is probed is still of interest. If 
the momentum factor is say 10, then even in going 
from 0.7 to 1.0 times the terminal speed, multi-line 
scattering is required of the radiation field seeping 
through the corresponding wind layers to account 
for the final acceleration. 

We have found that it is possible to obtain 
a simultaneous fit to both the continuum and 
the He II recombination lines from the ISO data 
for WR 136. We made calculations at two dif- 
ferent wind temperatures: Tw ~ 60000 K and 
Tw « 30000 K. Low clumping factors with Dc of 
1-3 are found. The wind velocity law is loosely 



constrained, but a value of 2-3 appears consistent 
with the various constraints imposed by the con- 
tinuum data, the emission line shapes, and the 
expectation of "turbulent" broadening in the flow 
arising from wind structure. Our estimate for (3 k, 
2 to 3 falls in the "gradual acceleration" range that 
is perhaps appropriate for multi-line scattering ac- 
celeration. 

For comparison a study of emission line pro- 
file variability by Lepine & Moffat (1999) included 
WR 136 as one of their sources. They observed 
He II 5411 and draw several interesting conclu- 
sions. First, they find that the; turbulence in the 
wind of WR 136 is typical of other WN stars in 
their study. And second, they derive the product 
» 75. For the radius of the wind base that 
we have adopted (i?* = 5Rq), this would imply a 
(3 of about 15. However, as the authors point out, 
as well as Koesterke, Hamann, & Urrutia (2001), 
exactly how one relates a /3- value derived from fea- 
tures observed in Doppler shift space to physical 
wind velocities is not clear, because there is almost 
certainly a flow of wind material past the clumps 
which are continually forming and dissipating in 
the wind. So although the data of Lepine & Mof- 
fat for WR 136 would imply an extensive region 
of wind acceleration, the f3 may not describe the 
ambient wind, nor does the i?* derived necessar- 
ily refer to the base of the wind. Our method is 
to infer /? through the influence of density on the 
formation of the IR continuum and He II lines. 
Our results are perhaps complementary to those 
of Lepine & Moffat, but we are not studying the 
same phenomenon. 

Our derived low level of clumping is also inter- 
esting, since it suggests that the large scale wind 
flow that is sampled by the radio observations and 
is also claimed to be unchimped (Nugis et al. 1998) 
extends fairly deep into the wind of WR 136 where 
the IR emission forms. Perhaps clumps that ex- 
ist deep in a wind can be dissipated by wind drag 
with the matter blending into the ambient wind. 

In future studies, one could attempt to use 
somewhat shorter near-IR spectral data to sample 
the wind acceleration at smaller radii; however, 
this will probably be limited by the fact that high 
electron optical depths and strong line-blanketing 
effectively "veil" the innermost wind acceleration 
zone. Perhaps there is no direct technique for ob- 
servationally probing the low velocity portion of 
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the densest WR winds. 
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A. Analytic Line Profile Shapes for Constant Expansion Velocity Winds 

The emergent intensity 1^ at impact parameter p within the Hne frequencies will be given by the sum of 
the line component 1^,1 and the continuum component I^^c- Adopting Sobolev theory for the line formation, 
assuming that the line optical depth scales with p^, and assuming the continuous emission and opacity is free- 
free, the respective intensities were given in eqs. (16) and (17). Implicit is that the wind is also isothermal. 
Further assuming that S^^i = B^,, the total intensity expression reduces to 

h{p) = 1^,1 + /,,e = B,[l- e-(-»-+-s)J ^ (Ai) 
which was also given in eq. (18). For rays that intercept the stellar core, the emergent stellar intensity is 

/.(p) = /.,,e-(^— +^«), (A2) 

where the Sobolev optical depth ts is taken to be zero for red-shifted frequencies in the line, since those 
resonance points all lie aft of the star. 

In the case of constant wind expansion, the line profile shape including the influence of the continuous 
opacity is analytic in a certain limit as we now show. First, the Sobolev optical depth is given by 



-4 -2 

TS = nx W 



n dw ,^ 9, W 



(A3) 



For constant expansion-velocity, this reduces to 



Ts = rix ^ (1 - /i^) = Ti sinep ^ (A4) 
where we have used the fact that p = xsinO. The free-free optical depth is given by 

= f (A5) 

For constant expansion- velocity, and re-expressing the integral in terms of p and 6, an analytic solution exists 
giving analytically to yield 

Tff = i Tc p"^ (6* - sin 6* cos 0) . ( A6) 

The value of Tmax in eqs. (Al) and (A2) depends on p as follows: 

_ /^^(e-sin^cos^) {p<l), 

Tmax — i TT Tc (p> 1) ^ ' 



2 ^ 



The observed flux of line emission at any frequency in the line will be 



= 27ri?,(T„)^ ^-^^Pdp [l-exp(-O.5r,M0)p-3)] 



+ / [l-exp(-0.5rc/i(^)p-3-rs sin^p'^)] 

+ J pdp [1 — exp(— O.Stttc p~^ — ts sin^p"^)] | , 



(A9) 
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where 



/i(6») =6'-sin6lcos6'. (AlO) 

For blue-shifted frequencies, Pmin = sin 9 since the isovelocity surfaces are cones of constant opening angle 
that intercept a spherical core photosphere. Although all of the stellar photosphere suffers absorption via 
the free- free opacity, only the annulus of the star from to 7r/2 will be diminished by the line opacity. At 
red-shifted frequencies, Pmin = 1 and the second integral in (A9) vanishes; there is no line emission for p < 1 
owing to stellar occupation. 

To solve these integrals, we make the change of variable = Tc/p^. The integrals become 



/. = 2nT^/^BAT^)^yj^^^_^ u-^du^^^[l-eM-0-5h{0)u')] (All) 



[ ° du j [1 - exp(-0.5/i(6l) - ts t"^ sin 6 u^)] 

j du [l - exp(-0.57r - ts sin 6u^)] > 



(A12) 



It is important to note that the frequency information is contained in the functions of 6 via = — /i, and 
also that these functions are constants of the integration which is over u. 

There is in fact an anal}d;ic solution to this form of integral given in Gradshteyn & Ryzhik (1994; Eq. 3.478, 
#2), for which the relevant result is 

poo . . 

y u-^ du (l - e-"" j = 1.5 r(4/3) a^/^. (A13) 

1/3 

This solution is applicable to the case at hand in the limit that Tc ^1, with the result 

/. = 3r(4/3) 7rB.(T^) T^/' + ^ yr^) . (A14) 
Interestingly, for a strong line with r; Tc, the profile shape reduces to 

U = 3r(4/3) 7rS,(T,„) rf/' (1 - wD'/^ (A15) 

The results of this derivation show that symmetric emission profiles result, which are overall similar in 
appearance to those observed for WR 136. The rounded shape of the profiles implies that a measurement 
of the widths of these lines would seriously underestimate the speed of the wind flow. Examples of the 
normalized profile shape of eq. (A15) for various values of the ratio ti/tc are given in Figure 4. 
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Fig. 1. — The geometry used in our calculations. 
All lengths are normalized to the stellar radius i?* . 
The observer is to the right. The observer coor- 
dinates are the cylindrical system (p, a, z), with 
azimuth a not shown. The star has spherical coor- 
dinates (r, d, (j)) . With the assumption of spherical 
symmetry, wc chose 6 measured from the z-axis 
as indicated. Also shown is a contour for an isove- 
locity zone with = constant and a ray passing 
through the envelope from the observer direction. 
At the frequency in the line corresponding to the 
Doppler shift Vz, the contribution to the observed 
emission along this ray includes line emission from 
the indicated point and continuum emission from 
both fore and aft of this point, along with the as- 
sociated attenuations. 
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Fig. 2. — Illustrates the influence of the velocity 
law on the optical depth (top) and emergent in- 
tensity (bottom) with Ii, = exp(— r,/,ff) along 
the line-of-sight to the star. The two cases shown 
are for /3 = (constant expansion-velocity) and 
/? = 1. Unlike a static exponential atmosphere, the 
extended envelope effects of an expanding wind 
lead to an attenuation profile that is not especially 
sharp as a function of radius. 
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Fig. 3. — Illustrates the influence of the veloc- 
ity law on the emergent line contribution function 
along the line-of-sight to the star plotted against 
the continuum optical depth (top) and wind ve- 
locity in the wind (bottom). At top, the two solid 
curves are optically thin lines for the indicated (3- 
laws. The dotted line is for an optically thick line 
with (3=1. The effect of large line optical depth 
is to reduce the amount of emergent line fltix and 
shift the peak location to larger radii. At bot- 
tom, just the (3=1 curves arc plotted against 
the normalized wind velocity. The range of veloci- 
ties over which line flux contributions arise is seen 
to be fairly broad. In this case, for the optically 
thin line, the peak in the line flux contribution oc- 
curs around v = 0.65i)oo and has a full width of 
w 0.24t;oo- 




Fig. 4. A sequence of theoretical emission pro- 
files for a constant expansion- velocity wind as the 
parameter t;/tc is varied, where n is the scale fac- 
tor for the line optical depth, and Tc is the scale 
factor for the free-free optical depth. The profiles 
are continuum normalized and plotted in velocity 
units normalized to the wind terminal speed. 




Fig. 5. The de-reddened low resolution 
(SWSOl) spectrum of WR 136. Major emission 
features, including the He II recombination lines 
used in the analysis in this paper are identified. 
Beyond about 6 fim, the signal-to-noise ratio is se- 
riously degraded, precluding derivation of the ve- 
locity law in the outermost regions of the stellar 
wind. 
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Fig. 6. — The de-reddened continuum flux distri- 
bution. This continuum is fit with a power law to 
give the index, = d log /a/^ log A = —2.87. 



1500 



„1000 

E 









v.= 1490 


km S-' 






T 

o 


o 






^ r- 


■ Hell III 


1 


1 1 




to 


■* 03 








- 




- 



t_J I I I I I 1 1 , \ 1 . I I I I I I -I 

2 3 4 5 6 



X(/Lim) 

Fig. 8. Velocity as a function of wavelength as 
determined by the HWHM of the He II recombina- 
tion lines. The terminal speed v^o = 1490 kms~^ 
is indicated as a dotted horizontal line. The dis- 
tribution is fit with a power law to give the index, 
7v = rflogHWHM/dlogA = 0.22. 
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Fig. 7. High resolution (SWS06) spectra of six 
He II recombination lines. The spectra have not 
been de-reddened. With reference to the red wing, 
the HWHM velocities measured for our analysis 
are indicated by vertical dashed lines. 




Fig. 9. — Shown as dotted lines are the Icr "al- 
lowed" zones from the ISO data. Model calcula- 
tions for the continuum slope (abscissa) and flux 
level (ordinate). Left is for a modest wind tem- 
perature of 27700 K and right for 57700 K. The 
three solid curves in each panel are for /? = 1, 
3, and 10. The points are for different clumping 
factors Dc = 1, 3, 10, 20, and 50, moving from bot- 
tom to top. The dashed lines are for models with 
a transition from He III to He II at r = 30i?*. 
Although the curves can shift significantly (espe- 
cially for lower temperature), the shift is mostly 
lateral such that low values of Dc are favored with 
or without the Helium ionization transition. 
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Fig. 10. — A sequence of model emission line 
profiles to illustrate the effects of varying dif- 
ferent parameters. At top, line profile calcula- 
tions with higher /3 have narrower half-widths. 
At middle, line profiles of higher optical depth 
Ti have stronger emission. At bottom, winds of 
higher temperature yield relatively stronger 
lines. (The same would apply for winds of 1ow(t 
Dc-values.) The lines are continuum normalized, 
so one should be aware that the variation of /3 and 
Tw changes the continuum. Evcni though these pa- 
rameters are not in fact independent, by varying 
them independently, we obtain a grid of models 
to illustrate the range of line widths and strengths 
that can result. 
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Fig. 11. — Application of the constant expansion- 
velocity line profile shape to the observed emission 
lines (plotted as continuum normalized). These 
theoretical asymptotic forms actually match the 
observed lines fairly well, although typically being 
somewhat broader (especially clear in the 14-11 
line). Even so, better fits arise when a velocity law 
is included, especially in light of the fact that there 
is some blending that likely adds extra broadening 
to the lines of interest. 
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Fig. 12. — A plot of the line broadening a re- 
quired to bring model line widths into agreement 
with that observed in the He II 7-6 line. The solid 
lines are for Tw = 57700 K and the dashed ones 
for Tw = 27700 K. In each case the clumping fac- 
tor Dc = \ for the upper curve and Z>c = 3 for 
the lower curve. The lower line is the minimum 
(J scit by the instrumental broadening. The two 
starred points indicate models with /3, Dc, and a 
that meet all the requirements fo the data, with 
the point at left for the higher temperature wind. 
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Fig. 13. — Shown as dashed lines are the pre- 
broadened line models with /? = 3, 15c = 2, and 
a wind temperature of Tw = 27700 K overplotted 
with the observed He II lines. 
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Fig. 14. — As in Fig. 13, but for models with 
[3 = 2, Dc = 2, and a wind temperature of 
Tw = 57700 K. 
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Table 1 
Journal of Observations 



Observation 


Spectrograph 


Date 


Observation 


Duration 


Name 






Number 


(seconds) 


WR 136 


SWS06 


1 Dec 1996 


38101710 


12872 


WR 136 


SWSOl 


1 Dec 1996 


38101711 


1140 


A\'R vm; oil 


SWSOl 


1 Dec 199(i 


.'-!<S101712 


1140 



Table 2: Line Identifications 



A Principal Line Blends at Av (km s 



Notes 



2.83 He II 9-7 C IV 15-9 at +750 



3.09 He II 7-6 



N III 8-7 at -390 

C III 9-8 at -290 
C IV 14-12 at 



features near the blue and red 
wings are possibly weak He I and 
He II lines 

the N line could be a significant 
contributor to the line flux 



He II 11-8 at +350 



since He II 14-9 have emission 
at about 30% of fcont, this feature 
likely has a slight influence on 
the width of He II 7-6 



3.15 He II 14-9 He II 19-10 at +560 



3.74 He II 16-10 H I (Pf 7) at 
He I 8-5 at 



4.65 He II 14-10 H I (Pf (3) at 

He II 22-12 at +1225 

4.76 He II 8-7 He II 17-11 at -1350 

7.21 He II 14-llW He II 19-13 at +1460 



blend may account for difficulty in 
matching the red wing width 



a sequence of transitions between 
different terms in the 8-5 line spread 
from —1000 to 0, which may account 
for the slight blue wing broadening 



"Wc include 14-11 so that there are 3 lines with upper level 14. The 14-9 and 14—11 do not blend with H, but 14—10 does. In 
comparing the two with odd lower level, we can confirm whether the lines are thin. By comparing to the one with even lower 
level, we can infer the H/He ratio. 
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Table 3: He ii Recombination Line Data for WR 136 



Wavelength 


Transition 


HWHM'^ 


Continuum Flux 


EW 


(/zm) 


(He II) 


(kms-i) 


(Jy) 


(A) 


2.826 


9-7 


1050 


3.44 


210 


3.092 


7-6 


1090 


3.20 


790 


3.145 


14-9 


1070 


3.17 


71 


3.739 


16 10 


1050 


2.79 


97 


4.652 


14-10 


1220 


2.37 


220 


4.764 


8-7 


1170 


2.32 


570 


7.208'' 


14 11 






180 



" The HWHM is evaluated for the red wing. 

^ The He II 14-11 line and neighboring continuum are not sufBciently reliable to quote a line width or 
continuum flux value, although it is shown for the model line fits described in §3.3.4. 



Table 4: Mass-Loss Rate Estimates for WR 136 



Reference 


M 


Waveband 


Flow Type 




(lO-^Meyr-i) 






Abbott et al. (1986) 


3.2 


Radio free-free 


smooth 


Barlow & Cohen (1975) 


2.3 


Radio free-free 


smooth 


Crowther & Smith (1996) 


12.5 


NIR lines 


smooth 


Hamann fc Koesterke (1998) 


12.5 


Optical/UV lines 


smooth 


Nugis et al. (1998) 


6.25 


Radio 


clumped" 



" These authors allow for clumping in their analysis, but conclude that the outer wind is unclumped. 
Table 5: Adopted Parameters for WR 136 



Parameter 


Value 


Subtype" 


WN6 




1.82 kpc 


rpc 


57700 K 




5i?0 




2.5 X lO^ig 


Md 


6.25 X IQ-^Mq yr" 




1490 kms-i 




27,700-57,700 K 


xf 


0.12 



" From van dcr Hucht (2001). 

^ From Lundstrom & Stcnholm (1984. 

" Ftom Crowther & Smith (1996) that are assumed for the base of the wind. 

From Nugis et al. (1998). 

From Paper I. 
^ X is the hydrogen mass-fraction. 
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